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ABSTRACT 

We present an analytical solution to the two-parabola Landau model, applied to melting of metal 
particles with sizes in the nanoscale range. The results provide an analytical understanding of the 
recently observed pseudo-crystalline phase of nanoscale Sn particles. Liquid skin formation as a 
precursor of melting is found to occur only for particles with radii, greater than an explicitly given 
critical radius. The size effect of the melting temperature and the latent heat has been calculated 
and quantitative agreement with experiments on Sn particles was found. 

Introduction 

Among the thermodynamical effects that depend on the size of the system, the depression of 
melting temperature when one or more dimension of a crystal is reduced to nanometer scale, is 
one of the best known. The melting point depression has been intensively investigated since it 
was experimentally discovered by Takagi HI in 1954. Melting of small metal particles (2j[3l|4j 
121 as well as semiconductor particles @ at temperatures below the bulk melting temperature is 
now experimentally well established. The phenomenon was already predicted theoretically by 
Pawlow in 1909. Since then several models have been proposed, for a review see |E1 15). 
Commonly these models predict a linear relationship between the depressed melting temperature 
and the inverse radius of the particle owing to the increase in the role of the surface free energy. 

More recently, studies of ultra-fine particles have indicated that small scale systems may ex- 
hibit quantum-size effects that modify the electronic states. Observations of morphology and in- 
ternal structures of ultra-fine particles show the presence of a pseudo-crystalline phase (where the 
particles are continuously fluctuating between different structures) BUS- The pseudo-crystalline 
phase has been discussed theoretically by ifTTI . 

The quantum-size effects occur when the number of surface atoms N s are comparable to the 
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number of bulk atoms N. For a spherical particle with radius R we have from lfT2l : 

N 2R 

where a is the atomic spacing of the crystal. A typical metal particle with radius R = 5nm will 
therefore have about 10% of the atoms at the surface. 

Premelting of surface layers has recently been studied intensively both theoretically and ex- 
perimentally. Applying Landau theory to a semi-infinite system, lfT3l [141 have shown that surface 
melting may induce critical phenomena. Later, |[T5l [T6l [T71 have studied surface melting of semi- 
infinite systems both experimentally and theoretically, and good agreement was found between 
experiments and the two-parabola Landau model. Surface melting of finite systems has been stud- 
ied experimentally by means of electron transmission microscopy (TEM) lfl2l [181 [191 and X-ray 
techniques dQlEll however only few theoretical studies have been reported ll22"Tl . 

To address the question of thermodynamical size effects and surface melting of finite size 
systems we take an analytic approach, solving the two-parabola Landau model suggested by lfI3Ti 
that was analyzed numerically for finite systems by 11231 . 

Model 

The Landau free energy functional F[M(R)] for a spherical particle of radius R described by the 
ordering parameter M(r) is given by 11231 

F[M(R)} = 4%R 2 f s (M) +4n£r 2 |/(M(r)) + ~ (^j j dr, (2) 

where the ordering parameter M = 1 describes the perfect crystal and M = describes pure liquid. 
The term f s {M(R)) is the free energy contribution from the surface per unit area and the term 
f(M(r)) is the free energy per unit volume for homogeneous bulk material. The second term of 
the integrand is the gradient in the order parameter and J is assumed to be constant [23 ]. The free 
energy function f(M(r)) is assumed to be polynomial and given by a two-parabola expression 

{%M 2 (r)+A(T) for M<M* 
2 (3) 
f(l-M(r)) 2 for M > M* 

where a is a material dependent parameter and M* is the intersection points of the two parabola: 

ir-I-*2>. ( 4) 

2 a 

A(r) is the difference in homogeneous bulk free energy per unit volume between the liquid and 
the solid phase at temperature T. Approximating to first order in the temperature 

A(T)^L h ^L (5) 
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where Ly is the bulk latent heat of melting per unit volume and T m is the bulk melting point, the 
surface energy term can be expressed as lfT4ll 

f s = -LM 2 (R)+ llv , (6) 

where a s is a material dependent parameter and y/ v is the liquid-vacuum interfacial energy per 
unit area. The parameters a s ,a and J are related to the interfacial energies YsZ>Ysv>Y/v an d Ay = 
Jsv — Jsl — Y/v, where y s j and y sv are the solid-liquid and solid- vacuum interfacial energies per unit 
area respectively. From lfT3Tl we have: 

I 

£, can be identified with the correlation length in the liquid state. The following relations are also 
given in lfl5Tl : 

4 1 + — 

oc = ^r J = *1s& a s = —^4y sl . (8) 

1 la 



Further K is defined as: 



Solution 



J 

%OL s 



(9) 



The order parameter M{r) that minimizes the total free energy F[M(R)] is found by applying the 
variational principle 8F[M}/&M = 0. The Euler-Lagrange condition then implies: 

d 2 M 2dM 1 . 

i^ + ;* + ? (1 - M) = a <10) 

The solution is given by a sum of the first and second kind of modified spherical Bessel func- 
tions of zero order: 

, , exp(— r/E) sinh(r/£) 
\-M{r)=A r ^ +B r K ^\ (11) 

where A and B are constants. The singularity at r = 0, implies A = while B is determined by the 
surface condition lfl4l : 

J d —{R) = ^h v (12) 
dr y ' dM(R) v ' 

The pure crystal phase then has an order parameter profile given by: 

1 r sinh(f ) 

M C RY(r) = l p : ttIt- (13) 



1+K 



coth(f)-| ^sinh(|) 



The crystal order parameter profile Mqry, equation [131 has been plotted for tin particles with 
different radii R in figure [T^. Naturally the atoms near the surface are less ordered than the interior 
of the particle. Surprisingly, significant disorder of the interior of the particle is seen for particles 
with radius R < 5^. We interpret this phenomenon to corresponds to the pseudo-crystalline phase 
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Fig ure 1: (a) Order parameter profile Mcry for Sn particles with different sizes. When R < 5t, a pseudo-crystalline 
phase occurs, (b) The order parameter profile for a Sn particle with radius R = lOi; as function of temperature. 

observed for ultra-fine tin particles by Oshima and Takayanagi @. The pseudo-crystalline phase 
has been observed for tin clusters with R < 2.5nm. Clusters where 2.5nm < R < 3.5nm were 
mainly observed to be crystalline while a small fraction were still in the pseudo-crystalline phase. 
The correlation length £ is typically of the order of Inm, thus giving good agreement between 
the model and the experiment of @. The origin of the pseudo-crystalline phase is, of course, the 
increasing importance of the number of surface atoms. 

For the sake of simplicity the order parameter profile M{r) including a quasi liquid layer (QLL) 
phase Mqll is given in the appendix. Figure £j})) shows the temperature-dependence of the order 
parameter profile M(r) for a Sn particle with radius R = 10^. The solid curves are for M(r) > M* 
and represent the solid state. The dashed curves are for M(r) < M* and represent the liquid and 
quasi liquid states. It is seen that the order profile M(r) does not change significantly due to the 
presence of the quasi liquid layer (QLL) and that the quasi liquid layer remains thin until the 
complete melting. 



Size dependence of the latent heat of fusion 

In classical thermodynamics the latent heat of fusion is assumed to be constant and independent 
of the size of the system. Recent experiments [24 J and MD simulations |[25ll26l . however indicate 
that the latent heat of fusion may be size-dependent. From the size effects of the order parameter 
profile M(r), see figure £[]), it is natural to expect a decrease in the latent heat of fusion as the 
particle size decreases. However it is not obvious how the order parameter profile can be linked to 
the latent heat of fusion. 
We suggest the following definition: 

, . 471 [*r 2 M 2 r]) y(r)dr 
L m {R)=L b Jo c fy (14) 
471 Jo r 2 Jr 



4 




20 30 40 

Radius R in nm 



Figure 2: Size dependence of the latent heat of fusion for Sn-particles, the datum-points are taken from 124] and the 
solid line is calculated from equation[T5] with % = \.5nm, y sv = 654mJ/m 2 , y; v = 570mJ/m 2 and y s / = 66mJ/m 2 . 



where Lb is the bulk latent heat of fusion. For a system with perfect crystal structure {Mqry (r) = 1) 
the latent heat of fusion L m (R) will be equal to the bulk latent of heat of fusion Ly. However the 
crystal structure of a finite size system is not always perfect, see figure (Q]). Surface effects makes 
the latent heat of fusion decrease with decreasing size because the order parameter M(r) is less 
than one as the particle surface is approached. With use of the crystal order parameter profile 
Mcry, equation (TOT ), the evaluation of the integral in equation (PT41 gives: 



L m {R)=L b {\ 



2 sinh 2 (7?/^) T 2 



1 | coth(/?A)[l-4r] 3^ | 6^ 2 



r 2 



2R FR 2 



where 



r= i+K[coth -z,/r]. 



(15) 



(16) 



Note that in the thermodynamical regime, R ^> the latent heat of fusion L m (R) is size inde- 
pendent, thus L m (R) = Lb. Throughout this paper the interfacial energies for Sn are taken as 
Ysv = 654mJ/m 2 , Y/v = 570mJ/m 2 and y s i = 66mJ/m 2 1 15 ] while the bulk latent heat of fusion for 
Sn is taken as Lb = 59.2 J/g ll27l or Lb = LbPe = 411 .2MJ/m 3 , where is the liquid density. The 
correlation length £ is estimated from ifTBI : 

where p/ = N /V is the particle number density of the liquid state and 1cb is the Boltzmann constant. 
For Sn p/ = 354 x 10 26 ra~ 3 and T m = 505K thus the correlation length ^ for Sn can be estimated to 
^ = \ .5nm. With these parameters equation ([TBI is compared with a experiment done on supported 
free Sn particles 11241 . in figure ©. Within the experimental uncertainty there is good agreement 
between the experimental and calculated data. 
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Figure 3: Size dependence of the melting temperature of supported free Sn-particles, the datum-points are taken 
from [24 1 and the solid line is calculated from equation [T9j with i; = \.5nm, y sv = 654mJ/m 2 , y/ v = 570mJ/m 2 and 
y,/ = 66mJ/m 2 . 

Sharp bulk melting 

It is known that ultra fine particles can not exhibit a strict phase transition. Finite-size scaling 
analysis shows that a first order phase transition of a finite system is broadened over an interval: 

(18) 

T c LN 

where kg is the Boltzmann constant, N is the number of particles and L is the latent heat per atom 
[28 ]. This broadened interval has also been found from quantum-statistical considerations by [29]. 
For a typical metal particle with a radius R = 5nm we find that AT C ~ O.IK, thus the transition is 
rather sharp. In figure \Vp, it is seen that the order parameter profile M(r)for a particle with dry 
surface does not change significant by surface melting. Therefore, it is a good approximation to 
consider melting of a small particle as a first order transition from a dry particle to a liquid particle. 
The melting transition occurs when the Landau free energy of the solid and the liquid phase equals. 
The critical melting temperature T c is therefore found by solving F[Mcry] = F[M = 0] where F is 
given in equation (©: 

T m -T c = 6y sl coth(fl/5)-5/fl 
T m RL b \ + K[coth{R/Q-^/RY 

In figure[3]the melting temperature T c of supported free Sn particles is compared with the model. 

With the same interfacial energy values as mentioned before good agreement was found except 

for the two smallest particle sizes. 

Remark that in the thermodynamical regime, R 3> ^, the melting temperature T c simplifies to the 
classical thermodynamical result: 

T m -T c _ 3(y sv -y lv ) 1 
T m Lf, R 
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where the depressed melting temperature (T m — T c )/T m is linear proportional with the inverse ra- 
dius R. In the opposite regime, R <C it was found from equation [l9]fhat the melting temperature 
T c is size independent: 



In this regime, a quantum mechanical model is necessary to describe the phenomenon properly. 
However the Landau model should be valid in the intermediated regime R « 

Surface and bulk melting 

Before the solid-liquid phase transition takes place, surface melting can occur as seen in figureQ}). 
When the intersection parameter M* exceeds the crystal order parameter Mqry surface melting 
starts. The temperature Tqll, for which the quasi liquid layer starts growing, is found by solving 
M* = McRy. Solving for Tqh to first order in R we get: 



For particles with large radius R the second term in equation d22l . caused by the curvature, will 
vanish and surface melting can only occur if Ay > 0. As the particle size decreases the curvature 
effect becomes more important. Generally, surface melting will occur when Tqll < T c and when 
the right hand side of equation (1221 is positive. In figure © a schematic phase diagram for surface 
melting is drawn. For Ay > there exists a critical radius R c , which is given by: 



The model predicts surface melting for particles with radius R larger than the critical radius R c , 
while small particles with radius R below the critical radius R c will maintain a dry surface. The 
existence of R c was also predicted from MD simulations on small Au particles by 11251 . They 
found that Au clusters containing less than 350 atoms, which corresponds to R c = Inm, do not 
show any quasi liquid layer. 

The critical radius R c has been measured experimentally for Pb 11301 and Sn PTl . The experiment 
of lf3Tl found R c /£, = 5.5 for Sn which is comparable with the model value R c /£, = 3.8. 

Surface premelting 

For particles with Ay > and R > R c , surface melting will occur before the transition from solid 
to liquid. By definition of the intersection point M*, the quasi liquid layer thickness lis found by 
solving M(r,l) = M* where M(r,£) is given in the appendix, equation d28l . The quasi liquid layer 




(21) 




(22) 




(23) 
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Surface non-melting 



Figure 4: Schematic phase diagram for melting of small particles. The region in which surface melting occurs is 
shaded. For definitions of symbols, see text. 



is found to grow with increasing temperature: 

S ln' 



t(T) 



1 



A(T)R 

2y s i 



(24) 



where a(R) is given in the appendix and A(T) is given in equation|5] It is seen that the quasi liquid 
layer thickness diverges when 



A(T) 



2jsi 
R 



(25) 



The model predicts however that the solid-liquid phase transition takes place before the quasi 
liquid layer diverges, see ngureQJ). Equation 25 is equivalent to the Gibbs-Thomson relation which 
in this model gives the upper temperature limit for which a solid particle can exist. Further the 
quantity 2y s j/R is known as the Laplace- or capillary pressure. The pressure of a particle is always 
greater than the surrounding pressure by the amount 2y s i/R. The quasi liquid layer therefore 
diverges since the difference between the particle pressure and surrounding pressure becomes to 
weak to keep the solid particle together. 



Discussion 

Surface pre-melting initiated from a fiat surface of semi-infinite systems studied by fT5l . is found 
in the limit R — > °°, thus equation (l24l simplifies to : 
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note that surface melting appears only for Ay > 0. For semi-infinite systems the quasi liquid layer 
thickness diverges for T — > T m . This is often interpreted as the bulk melting transition. However 
due to the finite size nature of any crystal the diverging quasi layer thickness is preempted by a 
first order melting transition. Take for example a film with thickness L, the melting temperature 
T c is given by P2l : 

A (r c ) = _L IkzJl = IflZll . ( 27) 

Thus the quasi liquid layer thickness at the melting temperature T c is of the order of £{T C ) ~ 
ln(LA). 

The present work has considered free finite size particles. However due to the generality of 
the model the obtained results could also be applied to coated particles or particles embedded in 
a matrix where both premelting |[T2l l33l and superheating P4l l35il have been observed. For such 
systems the surface energies y sv and Y/v should be replace by the interfacial energies, y m (solid- 
coat/matrix) and y,,, (liquid-coat/matrix). 

Conclusion 

In this contribution, an analytic solution to the two-parabola Landau model for finite size spherical 
systems was given. The solution provides insight into the size effect of the melting temperature 
and the latent heat of fusion. For particles with size comparable to the correlation length a non- 
linear dependence on size was found for both the melting temperature and the latent heat of fusion. 
For large systems compared to classical thermodynamic results was found. The model further 
provides an analysis of surface pre-melting. It was predicted that surface pre-melting only persists 
on systems with radius greater than a critical size. 
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Appendix 

The general order parameter profile M(r) including surface pre-melting is given in this appendix. 
For simplicity £, is taken as unit length. £ signify the quasi liquid layer thickness and £ is defined 
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by: I = R-l 



where 



M(r.l)-{ (28) 
M QLL {r) = ^ exp(r - i) {a + exp(2[/? - r])} for M < M* 



[coth(£)-i] [a + exp(2£)] 



exp(2£)(l + I)-a(l-i) 

(30) 



k(1 + I)-1 



and 



K(l-|)+1 



£ 1 

^~ l+£G + exp(2f)' (31) 
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